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1.  MICROMECHANICS  OF  DAMAGE 
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Micromechanics  consists  in  deriving  the  behaviour  of  materials  at  meso-scale  from 
the  study  of  specific  mechanisms  at  micro-scale.  The  micro-mechanisms  must  be 
precisely  defined  from  physical  observations  for  both  the  geometries  and  the  kinematics. 
Their  mechanical  modelling  is  performed  with  elementary  usual  constitutive  equations 
established  at  meso-  or  macro-scale  for  strain,  crack  growth  and  fracture.  When 
compared  to  the  direct  analysis  of  the  macroscopic  properties,  the  additional  power  of  this 
approach  comes  from  a  better  modelling  of  the  possible  interactions  between  different 
mechanisms  and  from  the  homogenization  that  bridges  the  gap  between  micro-  and  meso- 
scales. 

When  using  classical  continuum  thermodynamics  concepts,  the  effects  at  meso- 
scale  of  material  degradation  at  micro-scale  are  characterized  by  an  internal  variable  called 
damage.  Hereafter,  the  effects  of  the  material  degradation  are  assumed  to  be  isotropic  at 
meso-scale  and  the  damage  variable  to  be  a  scalar  denoted  by  D.  In  this  section, 
definitions  for  this  scalar  variable  are  derived  from  the  study  of  two  different  micro¬ 
mechanisms. 

1.1.  Brittle  isotropic  damage 

1.1.1.  Microcracks  and  scalar  damage  variable 

The  main  mechanism  of  brittle  damage  is  the  nucleation,  growth  and  coalescence 
of  microcracks  up  to  the  initiation  at  meso-scale  of  a  crack.  Hereafter,  a  relationship 
between  the  micro-crack  pattern  and  the  damage  variable  D  is  established. 

Let  us  consider  a  Representative  Volume  Element  at  meso-scale  as  a  cube  of 
dimension  (1  *  1  *  1).  This  RVE  is  assumed  to  be  constituted  at  microscale  of  cubic  cells 
of  dimension  (d  *  d  *  d)  in  which  may  lie  a  microcrack  of  any  area  Si  and  any  orientation 
(see  Fig.  1).  The  number  of  cells  is  m=P/d3  and  the  number  of  cracks  n  <  m. 

/  / 


Fig.  1.  Micro-  and  meso-models  for  brittle  damage 
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The  geometry  being  defined,  the  modelling  consists  in  writing  the  balance  of  the 
dissipated  energy  calculated  by  classical  fracture  mechanics  one  one  hand,  and  calculated 
by  continuum  damage  mechanics  on  the  other  hand. 

For  a  cracked  cell  i,  subjected  to  a  given  state  of  stress,  if  Gi  denotes  the  strain 
energy  release  rate  corresponding  to  a  crack  of  area  sj,  Di  the  equivalent  damage  of  the 
cell  and  Yj  the  strain  energy  density  release  rate,  the  balance  of  dissipated  energy  can  be 
written  as  : 

G,  sj=  Yi  Di  d3 

For  the  n  cracked  cells  of  the  meso-cube,  the  previous  relation  becomes  : 
n  .  n 

I  Gi  Si  =  2  Yi  Di  d3 

1  1 

Assuming  that  brittle  growth  of  microcracks  occurs  at  G  =  Gc  =  constant,  corresponding 

to  Y  =  Yc  =  constant,  it  can  be  written  that : 

n  .  n  . 

Gc  I  Si  =  Yc  d3  I  Di 
1  1 

Furthermore,  if  it  is  assumed  that  Sj  =  0  corresponds  to  D;  =  0,  integration  of  the  previous 
relation  yields : 

Gc  I  Si  =  Yc  d3  I  Di 
1  1 


The  simplest  homogenization  consists  in  defining  the  damage  D  at  macro-scale  by 

the  mean  value  of  the  damages  D,  of  all  the  l-'/d3  micro-cells,  so  that 
1  m  !  n  d3  n 

D=^D‘  =  ^D‘  =  wI  D‘ 


or 


1 

r,„  n 

Si 


l 


i 


The  term 


Gc 


"r  may  be  evaluated  from  a  rupture  criterion  defining  the  initiation 
n  Y  c  1- 

of  a  meso-crack.  Because  of  the  localization  of  the  damage  phenomenon,  it  can  be 
assumed  that  the  meso-crack  initiation  occurs  when  only  part  of  the  flat  volume  (1  *  1  *  d) 
is  micro-cracked,  the  other  microcraks  being  neglected. 

In  other  words,  if  it  is  assumed  that  the  mesocrack  initiates,  i.e.  D  =  Dc  at  meso-scale, 
when 


then, 


I  sj  =  k  12 
1 


dc  =  TT^3  kl2 


Gc 

Ycl3 


kl2 
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Hence, 


2  s, 

1 


Dc 

k 


n 

In  this  case,  the  damage  variable  D  appears  as  the  micro-cracks  surface  density  (X  s,  /  l2) 

1 

corrected  by  a  factor  (here  Dc/k). 

If  the  following  simplest  fracture  criterion  is  considered 

n 

I  Si 

x  Sj  =  l2  —>  k  =  1  — »  Dc  =  1  then,  D  = 

l 


By  the  way,  this  calculation  gives  an  order  of  magnitude  of  a  characteristic  length 
which  allows  for  the  matching  between  fracture  mechanics  and  damage  mechanics  namely 
1,  the  size  of  the  Representation  Volume  Element.  Since, 

D  _  Gc  _  Dc 
"  =  Yc  1*  ~  k  l2 

I  Si 
1 

1  =  k  Dc  Vs 

1  0 

Hence,  for  the  simple  fracture  criterion  k  =  1,  Dc  =  1, 


For  most  metallic  materials 


:ht  alloys 

Steel  and  high  alloys 

0.005 

<GC< 

0.05 

MPa  m 

2. 

<YC< 

10. 

MPa 

0.0025 

<  1  < 

0.005 

m 

and  for  concrete  in  tension  Gc=  3  10^  Mpa.m,  Yc  =  1.5.10-4  MPa,  so  that  1  s  2.10'E 
This  shows  that  the  size  of  the  physical  RVE  must  be  of  the  order  of  the  millimeter  for 
metals  and  of  the  order  of  the  decimeter  for  concrete. 
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1.1.2.  Brittle  damage  growth 


As  a  specific  example,  let  us  derive  the  kinetic  damage  evolution  law  at  meso-scale 
which  corresponds  to  the  fatigue  microcracks  growth  at  micro-scale  of  Fig.  1.  For 
simplicity  sake,  the  analysis  is  restricted  to  a  two  dimensional  problem,  for  which  e 
denotes  the  uniform  thickness  of  the  whole  RVE. 


With  Dc  =  k,  it  has  been  established  that : 
n  n  . 

S  Si  I  Si 

D=-W-  and  D-1 
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and,  for  each  cracked  cell,  if  2a  denotes  the  crack  length,  the  following  relations  hold  : 
Ki  =  (EG;)17"  Kj  =  ^  Gf1/2  G{  sv  =  2e  a, 


The  surface  growth  rate  s;  of  each  crack  can  be  expressed  as  a  function  of  the  strain 
energy  release  rate  Gj  of  the  corresponding  cell,  by  means  of  the  Paris'  law  of  fatigue 
crack  growth.  If  N  denotes  the  number  of  cycles  of  loading  in  mode  I,  and 
K\t  =  KMax  the  amplitude  of  the  stress  intensity  factor  (with  Kmin  =  0),  then 


5a 

5N 


=  C  K71 


w  here  C  and  r]  are  two  material  constants,  with  rj  h  4  for  many  metallic  materials. 

If  one  assumes  that  this  Paris'  law  corresponds  to  the  integration  over  one  cycle  of : 

a  =  r)  C  KT'  K  , 

then 

•  ri/o  • 

si=  CeE^Gj2  Gi 


A  relationship  between  Gj  and  Yj  can  be  found  through  their  definition  from  the 
elastic  energy.  If  w  denotes  the  elastic  strain  energy  density  and  Wj  the  elastic  strain 
energy  of  the  elementary  cell,  then  : 


G, 


whereas 


Yi  = 


3wi 

5d[ 
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Since 


In  this  example,  the  damage  rate  is  an  increasing  function  of  the  strain  energy  density 
release  rate  ;  for  most  materials  since  rj  =  4,  the  damage  rate  is  quasi-proportional  to  the 
strain  energy  density  release  rate.  The  damage  rate  is  also  proportional  to  the  rate  Y. 
This  will  be  used  in  section  2.1  as  a  guideline  to  derive  a  general  kinetic  law  for  damage 
evolution. 
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1.2.  Ductile  isotropic  damage 

1.2.1.  Microcavities  and  scalar  damage  variable 

The  main  mechanism  of  ductile  damage  is  the  nucleation,  the  growth  and  the 
coalescence  of  microcavities  by  large  local  plastic  deformations.  Hereafter,  a  relationship 
between  the  density  of  micro  voids  and  the  damage  variable  D  is  established. 

Let  us  consider  again  a  Representative  Volume  Element  at  meso-scale  as  a  cube  of 
dimension  (1  *  1  *  1).  This  RVE  is  assumed  to  be  constituted  at  microscale  of  cubic  cells 

of  dimension  (d  *  d  *  d)  in  which  may  lie  a  void  of  volume  d3  (Fig.2.). 

/  / 

/  / 


On  this  very  simple  geometry,  the  modelling  consists  in  writing  the  balance  of  the 
dissipated  energy  calculated  from  the  growth  of  the  cavities  on  one  hand,  and  calculated 
by  continuum  damage  mechanics  on  the  other  hand. 

For  the  geometrical  model  under  consideration,  the  porosity  P  can  be  defined  as  : 


where  p  and  po  are  the  current  and  initial  porosity  respectively,  and  n  the  number  of 
cavities. 

According  to  Gurson's  model,  the  porosity  P  at  meso-scale  is  equal  to  the 
hydrostatic  part  of  the  plastic  strain  due  to  the  growth  of  voids. 

Here,  this  assumption  leads  to  the  following  equality  written  at  meso-scale  : 
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•p 

At  meso-scale,  given  an  homogenized  stress  Cjj  and  a  plastic  strain  rate  e^,  'he  total 
power  dissipated  is 

This  can  be  split  in  two  parts  by  means  of  the  deviatoric  and  hydrostatic  quandries  i.e. 

^  =  (aijD  +  aH5ij)(efjD+  ePH  S(J) 

D  *pD  _  *p 

=  a  e  +  3  au  eu  +  0 

ij  ij  H  H 

The  first  term  is  the  power  dissipated  in  pure  plasticity  by  slips.  The  second  term  which 
corresponds  to  the  irreversible  change  of  volume  may  be  interpreted  as  the  power 
dissipated  in  the  RVE  to  increase  the  material  discontinuities  by  growth  of  the  mvities. 
This  latter  part  must  balance  the  damage  dissiparion  i.e.  : 


3  aH  1?  =  Y  D  l3 


so  that 


D  = 


3  a, 


3  oL 


LI 

Assuming  for  simplicity,  proportional  loading,  perfect  plasticity,- -y—  =  const,  and  the 

initial  condition  P  =  0  D  =  0.  the  integration  yields  : 

3  cLt  (p 


D  =- 


'H 


Y  “  p 

As  for  brittle  damage,  because  of  the  localization  of  the  damage  phenomenon,  it  can  be 
assumed  that  the  meso-crack.  initiation  occurs  when  a  set  of  cavities  occupies  only  pan  of 
the  flat  volume  (1  *  1  *  d),  the  other  cavities  in  the  RVE  being  neglected. 


In  other  words,  the  critical  value  of  the  porosity  corresponding  to  D  =  1  is  assumed  as 

d 
1 


p  _  n  w  lid  _  d 
rc  -  n  p  -  K  p  -  K 


3  o 


H 

This  allows  for  the  calculation  of  the  term— y~"  in  the  damage  equation  : 
3  j 


D  =  1  = 


Y 


-  k  t 


D  = 


n  d2 


1 1 


Hence, 


1.2.2.  Ductile  damage  growth 


As  a  specific  example,  let  us  derive  the  kinetic  damage  evolution  law  at  meso-scale 
which  corresponds  to  the  voids  growth  at  micro-scale.  For  simplicity  sake,  the  analysis 
is  restricted  to  the  particular  case  k  =  1 . 

The  kinetic  law  for  damage  evolution  can  be  directly  derived  from  the  expression 
for  D  established  in  the  previous  section,  so  that 


•  d-  • 

D  =  pr  n  +  2  n 


d  d 
P 


The  first  term  accounts  for  the  increase  of  the  number  of  cavities  and  n  denotes  the 
number  of  cavities  nucleated  per  second.  The  second  term  accounts  for  the  cavity 
growth.  In  the  Gurson  model,  the  porosity  rate  is  also  the  sum  of  two  terms  accounting 


for  nucleation  and  growth. 


a)  Damage  growth  by  nucleation  of  cavities 


To  model  nucleation.  Tvergaard  proposed  the  following  kinetic  law  for  porosity 

1*  =  A  Cl, ,  +  B  a,, 

eq  H 

u  here  A  anti  B  are  material  parameters. 


Assuming  for  simplicity  sake  a  sudden  nucleation  of  cavities  of  a  fixed  size  d 


and 

••11  •  • 

D  -  P  7  =  7  (A  o  +  B  o  ') 
d  d  eq  H 

so  that. 


D  =  d 


O I 


(A  +  B  t*1) 

C5oq 


Damage  can  be  expressed  as  a  function  of  the  accumulated  plastic  strain  rate  p  written  in 
terms  of  the  plastic  tangent  modulus  Et.  Assuming  proportional  loading,  i.e. 


oh 

OH 

• 

— 

Ooq 

°oq 

• 

2 

*P 

P  = 

3 

eu 

£  =  c 

ij  Ey 


so  that. 


D=  jEt 


'a  +  b 

_ _ Opq 


12 


b)  Damage  growth  by  enlargement  of  a  fixed  number  n  of  cavities. 


The  problem  of  void  growth  has  received  much  attention  in  the  past  20  years. 
Essential  results  are  the  McClintock  and  Rice  &  Tracey  analyses  which  derive  the  rate  of 
growth  of  a  cylindrical  or  spherical  cavity  of  volume  V  in  a  perfectly  plastic  infinite  body 
as  a  function  of  the  accumulated  plastic  strain  rate  p  and  triaxiality  ratio  On/Ccq  ,  i.e. 

♦  •  'X  CTij 

V  =  0.85  V  p  exp(  j  —  ) 

z  Obq 

Taking 

V  =  d-' 


leads  to 


•  •  ^  CTl-? 

3  d2  d  =  0.85  dJ  p  exp(  f  — ) 


-  a, 


Since 


dd  „  d2  ‘  d 
D  =  2  n  -p-  ,  D  =  n  ,  D  =  2Dd 


D  =  0.57  D  p  expl 


a 


°d) . 


c)  Conclusions 

In  this  example,  the  damage  rate  by  nucleation  and  growth  of  voids  appears  as 

-  proportional  to  the  accumulated  plastic  strain  rate, 

-  an  increasing  function  of  the  triaxiality  ratio  — , 

°cq 

-  through  Ej  or  D,  a  function  of  the  current  state  of  the  material. 

This  will  be  used  in  section  2.1  as  a  guideline  to  derive  a  general  kinetic  law  for  damage 
evolution. 
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GENERAL  PROPERTIES  AND  FORMULATION 
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2.1.  A  general  form  for  continuum  damage  evolution  law 


a)  General  formalism 


Within  the  framework  of  classical  continuum  thermodynamics,  the  thermo¬ 
mechanical  state  of  a  material  is  described  by  the  following  set  of  independent  state 
variables : 

V  =  (£,  T,  V)  with  V  =  (D,  £p,  Vp), 
where  £,  £p  denote  the  total  and  plastic  strain  tensors  respectively,  T  the  temperature,  V 
the  set  of  the  internal  variables,  and  D  the  damage  variable.  For  simplicity  sake,  hereafter 
only  isothermal  situations  will  be  considered. 


The  reversible  behaviour  is  described  by  the  Helmholtz  specific  free  energy 
H'  =  4/(£,  V), 
chosen  as 

vp  =  vpc(S  _  SP5  d)  +  4,P(SP,  Vp) 

with 

\  (S-  £P):(1—  D)30:(S  —  £P) 

where  denotes  the  elasticity  matrix  of  the  undamaged  material.  The  thermodynamic 
forces 


A  =  Ur,  A) 
'om  t 
3VF 


are  defined  from  the  following  state  laws  : 


3-; 


A  h  (-  Y,  -  c,  Ap)  =  p 


3^ 

av 


where  p  denotes  the  mass  density,  Z  the  stress  tensor,  and  Y  the  strain  energy  density 
release  rate  or  damage  energy  release  rate. 

The  irreversible  behaviour  is  described  by  a  dissipation  potential 

0>  -  0(3,  A  ;  V)  , 


from  which  the  following  evolution  laws  are  derived  : 


£P  =  - 


30 


V  =  — 


30 

3A 


In  particular,  the  damage  evolution  law  can  be  chosen  such  that 


b)  damage  vs  (micro -)plasticity 


In  section  I.,  it  has  been  established  that  damage  is  always  related  to  some 
irreversible  strain  either  at  micro-  or  meso-level.  This  property  can  be  taken  into  account 
in  the  evolution  law  for  the  damage  variable  by  assuming  that  the  factor  X  is  proportional 
to  the  accumulated  plastic  strain  so  that 


D  = 


5Fp 

9Y 


P 


The  irreversible  nature  of  damage  is  directly  taken  into  account  by  the  fact  that  the  variable 
p  is  always  positive  or  null. 


In  most  materials,  a  certain  amount  of  plasticity  must  be  accumulated  before 
damage  at  meso-level  appears.  In  metals,  this  corresponds  to  the  accumulation  of  micro- 
stresses  in  the  vicinity  of  initial  defects,  of  dislocations,  ...,  prior  to  the  nucleation  of 
micro-cracks  or  micro-voids.  To  model  this  phenomenon,  since  the  damage  evolution  is 
governed  by  the  accumulated  plastic  strain  rate,  it  is  natural  to  introduce  a  threshold  pp  on 
the  variable  p,  such  that 

D  =  §  P  if  P  -  Pd 

D  =  0  if  p  <  po 

or 

D  =  ^  p  H(p-pD) 
where  H()  is  the  Heavyside  step  function. 

In  monotonic  loading,  pp  can  be  identified  as  the  uniaxial  damage  threshold  e^,  whereas 
for  fatigue  or  creep  processes  pp  is  a  function  of  the  applied  stress,  as  it  will  be  discussed 
in  section  2.2. 


c)  driving  force  for  damage 

From  the  thermodynamical  analysis,  it  has  been  deduced  that  the  driving  force  for 
damage  is  the  strain  energy  density  release  rate  Y.  Hence,  Fp  must  be  a  function  of  Y  : 

Fp  =  FD  (Y,  ...) 
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d)  influence  of  the  triaxiality  ratio 


Another  important  feature  of  fracture  mechanisms  is  the  influence  of  the  triaxiality 
ratio  (CTH/Ceq).  where  it  is  recalled  that  oh  denotes  the  hydrostatic  stress  and  Ocq  the  Von 
Mises’  equivalent  stress.  This  effect  is  directly  taken  into  account  through  the  damage 
energy  release  rate  Y  which  is  a  function  of  the  triaxiality  factor  Rv  : 

Y  _  Peg2  Rv 
2  E  (1  -  D)2 

with 

O  /  (J  LT  \2 

Rv=  T  (1+V)  +  3(1-  2v)  -s- 

^  l  °eq  )  . 


e) a  general  form 

In  order  to  choose  the  proper  and  simplest  expression  for  Fq,  let  us  recall  the 
kinetic  damage  laws  obtained  by  micromechanics  for  particular  mechanisms  in  section  1.  : 


Brittle  damage  by  fatigue  growth  of  micro-cracks  : 


D  T1  C  Eh/2  d7172  e 


r|  being  of  the  order  of  4 

D  =  (const)  *  Y  Y 

Although  no  plasticity  has  been  introduced  in  the  analysis,  it  always  exists  at 
micro-scale  at  the  crack  tips  of  the  micro-cracks  and  it  is  possible,  at  least  formally 
to  relate  Y  to  p  through  a  plasticity  constitutive  equation  : 

Y  =  Y(<JCq)  -»  Y(acq) 
and 

•  • 

acq  (p)  ->  acq(p) 


Ductile  damage  by  nucleation  of  micro-cavities  : 


( 

A  +  B 
V 


Oh" 


P 


Ductile  damage  by  enlargement  of  micro-cavities  : 


D  =  0.57  D  p  exp  | 


ojA 

<w 
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The  qualitative  conclusion  which  can  be  drawn  from  these  three  results  is  that  the 
damage  rate  D  can  be  considered  as  proportional  to  Y,  which  is  a  function  of  (CH/CTeq), 
and  p . 

D-Yp 

or 

Fd  ~  Y- 

As  in  any  realistic  constitutive  equation,  a  material  dependent  scale  factor,  such  as 
(const),  ^  Ex,  or  0.57,  must  be  introduced.  Let  us  denote  by  S  this  material  constant  so 

that 


Finally,  according  to  the  qualitative  properties  listed  above,  the  damage  potential  is 
naturally  written  as : 


Y- 

Fd(Y;  (p.D))  =  xg  Hl'p-po) 

where  the  factor  2  has  been  introduced  to  compensate  for  the  factor  (1/2)  coming  from  the 
derivation.  Hence,  the  proposed  general  continuum  damage  evolution  law  is  the 
following  : 


I  D  =  |  p  H(p-p„) 


where  two  material  dependent  parameters  are  introduced,  viz.  S  and  pD  which  characterize 
the  energetic  resistance  against  the  damage  process  and  the  damage  threshold, 
respectively.  The  effects  of  the  temperature  T  are  taken  into  account  through  the  variation 
of  these  coefficients  with  T  and  through  the  accumulated  plastic  strain  rate  p  which  is  also 
a  function  of  T. 

Several  important  properties,  though  not  directly  introduced  in  the  formulation, 
are  also  naturally  exhibited  by  this  general  evolution  law,  i.e.  : 

-  the  non  linear  accumulation  of  damage, 

-  the  effect  of  mean  stress  in  fatigue, 

-  the  non  linear  interaction  of  different  kinds  of  damage. 
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2.2.  Damage  threshold 


Under  monotonic  loading,  the  damage  threshold  po  can  be  identified  with  the 
uniaxial  damage  threshold  whereas  in  the  case  of  fatigue  or  creep  loadings  it  is  a 

function  of  the  applied  stress.  It  corresponds  to  the  critical  level  of  plasticity  which 
induces  the  nucleation  of  microcracks  without  any  consequence  on  the  mechanical 
properties  and  can  be  related  to  the  energy  stored  in  the  material. 

Experiments  in  fatigue  have  shown  that  the  total  plastic  strain  energy  dissipated 

may  reach  tremendous  values  before  failure  although  the  stored  energy  remains  constant 

at  microcrack  initiation.  This  stored  energy  is  the  result  of  microstress  concentrations 

which  develop  in  the  neighbouring  of  dislocation  networks  in  metals  and  of 

inhomogeneities  in  other  materials.  For  a  unit  volume,  it  is  equal  to  the  difference 

t 

between  the  total  plastic  strain  energy  (  j  Cjj  e,j  dt )  and  the  energy  dissipated  in  heat  as 

0 

given  by  the  Clausius-Duhem  inequality  of  the  second  principle  of  thermodynamics. 


For  instance,  in  the  case  of  a  material  exhibiting  kinematic  X  and  isotropic  R 
strain  hardenings  and  no  damage,  under  an  isothermal  transformation  the  rate  of  energy 
dissipated  in  heat  is  : 

0  =  oij  efj  -  R  p  -  Xy  ctjj  >  0 

This  expression  may  be  calculated  from  : 

-  the  potential  of  dissipation  e.g. 

F=  (tr0  -  X)Cq  —  R  —  Gy  —  Xy  Xjj 


-  its  associated  normality  flow  rule 

•n  r)F  •  •  •  • 


da  i 


3R 


Ctii  =- 


_5F_ 

9X;j 


-  and  the  yield  criterion 
f=  (cD-X)eq  -  R  -  cy  =0 


so  that 

0  =  (Qy +  2^  XU  X ij)P 

Hence,  the  stored  energy  Ws  as  a  function  of  time  t  is  : 


t  »n  t  /  3  \  • 

Ws(t)=  |  Cjj  ejj  dt  -  r  CTy  +  j~x~  X'j  x  i j ) P  dt 
0  ok  °°  J 
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This  formula  can  be  simplified  if  the  following  assumptions  are  made  : 

-  the  effect  of  the  kinematic  hardening  is  neglected  : 

Ws=  j  (aeq  -  ay)  dp 

o 

-  the  variation  of  <Teq  is  neglected  as  for  a  quasi  perfectly  plastic  material 
Ws  =  [Sup  (<Jeq)  “  CTy]  P 

If  this  stored  energy  is  considered  as  a  constant  for  the  damage  threshold  po,  its 
value  can  be  identified  in  the  one-dimensional  monotonic  case  used  as  a  reference  with 
Pd  =  ePD- 


In  the  particular  case  considered  above,  and  with  the  crude  approximations  made, 
the  onset  for  the  damage  process,  or  damage  threshold,  corresponds  in  the  case  of  a 
monotonic  loading  to  the  ultimate  stress  ctu  : 

[Sup  —  <Jy]  pO  =  (0U  —  Oy )  Epjy 

so  that 

P  d  —  epn 

_ l)  Sup(a0M)-gv 


As  a  summary,  in  the  particular  case  considered  in  this  section,  the  whole  set  of 
equations  that  governs  the  damage  evolution  is  : 


A  V  •  It,  _  ou-ov 

D  =  c  P  H(  p  -  ep  - — - — - ) 

5  Sup(oct,)-ov 
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3.  LOCAL  FAILURE  CRITERIA 
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The  ultimate  stage  of  continuum  damage  evolution  corresponds  to  meso-crack 
initiation.  As  a  first  approximation,  this  critical  damage  state  can  be  characterized  by  a 
critical  value  Dc  of  the  damage  variable,  for  instance  Dc  =  1. 

In  fact,  the  value  of  Dc  is  not  solely  material  dependent.  As  discussed  below,  it 
also  depends  on  the  local  stress  state  . 

It  is  proposed  to  relate  the  conditions  for  the  meso-crack  initiation  to  the 
conditions  for  the  localization  of  the  deformation  in  the  material.  In  this  section,  these 
latter  conditions  are  studied  in  details,  inside,  at  the  boundaries  or  at  the  interfaces  of  rate- 
independent  solids  for  both  the  linear  and  the  non-linear  case.  Physical  interpretations  of 
these  conditions  are  also  given. 

3.1.  Introduction 

Since  the  pioneering  works  of  Hadamard,  Hill,  Mandel  and  Rice,  the  localization 
of  the  deformation  in  rate-independent  materials  is  treated  as  the  bifurcation  of  the  rate 
problem. 

Here,  we  give  a  general  view  of  the  various  bifurcation  and  localization 
phenomena  for  possibly  heterogeneous  solids  made  of  rate-independent  materials.  Under 
the  small  strain  assumption,  the  behaviour  of  these  materials  is  described  by  the  following 
piece-wise  linear  rate  constitutive  laws  : 

L  =  3  when  /  <  0,  or  /  =  0  and  b  :  S  :  e(v)  <  0 

( 1 )  c  =  1 :  E(v)  with 

L  =  K  when  /  =  0  and  b  :  5  :  £(v)  >  0 

•  • 

where  C  and  S(v)  =  s  respectively  denote  the  stress  and  strain  rates,  v  the  velocity,  and/ 
the  yield  function. 

The  general  class  of  materials  modelled  by  relations  (1)  includes  for  instance  the 
elasto-plastic  damageable  solids  the  behaviour  of  which  is  described  by  the  constitutive 
equations  discussed  in  section  2.1. a. 
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It  is  shown  that,  in  general,  different  types  of  localization  phenomena  may  occur, 
depending  on  the  failure  of  one  of  the  three  conditions  which  are  described  in  section  3.2. 
Their  physical  interpretation  is  the  following  : 

-  the  ellipticity  condition  is  very  classical.  Its  failure  is  the  condition  for 
localization  given  by  Rice  and  linked  to  the  appearance  of  deformation  modes  involving 
discontinuities  of  the  velocity  gradient.  It  has  also  been  related  to  stationary  acceleration 
waves  ; 

-  the  boundary  complementing  condition  governs  instabilities  at  the  boundary  of 
the  solid.  Its  failure  leads  to  deformation  modes  localized  at  the  boundary  and  is  related 
to  stationary  surface  waves  (for  instance  Rayleigh  waves) ; 

-  the  interfacial  complementing  condition  governs  instabilities  at  interfaces.  Its 
failure  leads  to  deformation  modes  localized  at  each  side  of  the  interface  and  is  related  to 
stationary  interfacial  waves  (Stonelv  waves). 

3.2.  Rate  problem  analysis  :  the  linear  case 

Let  us  consider  for  instance  the  body  sketched  in  Fig. 3. 


Fig.  3.  Different  types  of  localization  modes 
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Qualitative  results  can  be  exhibited  from  the  analysis  of  the  rate  problem  for  the 
so-called  linear  comparison  solid  (see  Hill).  In  this  case,  this  linear  problem  is  well- 
posed  if  and  only  if  the  following  conditions  are  met : 

-  the  ellipticity  condition  :  the  rate  equilibrium  equations  must  be  elliptic  in  the 
closure  of  the  body  Cl,  i.e. 

det(n  .  >1 .  n)  *  0  for  any  vector  n  *  0,  and  any  point  M  e  Q. 

-  the  boundary  complementing  condition  :  this  relation  between  the  coefficients  of 
the  field  and  boundary  operators  must  be  satisfied  at  every  point  P  belonging  to  the 
boundary  T  where  the  boundary  conditions  are  formally  written  as  3(v)  =  g.  This 
condition  is  easily  phrased  in  terms  of  an  associated  problem  on  a  half  space  defined  by 
z  >  0.  It  requires  for  every  vector  k  =  (kj,  k2,  0)  *  0,  that  the  only  solution  to  the  rate 
equilibrium  equations  with  constant  coefficients  (equal  to  those  of  the  operator  at  point 
P),  in  the  form 

v(x,  y,  z)  =  \v(z)  expfi  (kj  x  +  k2  >')] 

with  bounded  \v  and  satisfying  the  homogeneous  boundary  conditions  3(v)  =  0.  is  the 
identically  zero  solution  v  =  0. 

-  the  interfacial  complementing  condition  :  this  relation  between  the  coefficients  of 
the  field  operators  in  Qi  and  Q2  must  be  satisfied  at  every  point  Q  of  the  interface  1 
between  Qi  and  fb-  This  condition  is  again  easily  phrased  in  terms  of  an  associated 
problem  on  the  whole  space  divided  by  the  plane  interface  z  =  0.  It  requires  for  every 
vector  k  =  (kj.  k2,  0)  *  0,  that  the  only  solution  to  the  rate  equilibrium  equations  with 
constant  coefficients  (equal  to  those  of  the  operators  at  point  Q.  in  Q]  for  z  <  0  and  in  fb 
for  z  >  0),  in  the  form 

( v  1  ( x ,  y,  z),  v2(x,  y,  z))  =  (wi(z),  W2(z))  exp[i  (kj  x  +  k2  y)] 

with  bounded  (>vi,\vf2)  and  satisfying  the  continuity  requirements  (continuity  of  the 
velocity  and  the  traction  rates)  across  the  interface  z  =  0,  is  the  identically  zero  solution 
(\q(z),  v'2(z))  s  (0,  0) ;  (where  vq  and  vy  are  the  solutions,  respectively  for  z  <  0  and  for 
z  >  0). 


When  these  three  conditions  are  fulfilled,  the  rate  boundary  problem  admits  a 
finite  number  of  linearly  independent  solutions,  which  depend  continuously  on  the  data, 
and  which  constitute  diffuse  modes  of  deformation. 
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Remarks  : 

-  these  three  conditions  are  local,  and  this  is  particulally  important  when 
considering  their  numerical  implementation  ; 

-  the  above-given  results  remain  valid  for  an  arbitrary  number  of  non-intersecting 
interfaces,  an  interfacial  cordition  being  written  for  each  interface  ; 

-  the  failure  of  these  conditions  can  be  interpreted  as  localization  criteria  as  recalled 
in  section  3.1.  These  localization  criteria  can  also  be  used  as  indicators  of  the  local  failure 
of  the  material ; 

-  both  boundary  and  interfacial  complementing  conditions  fail  in  the  elliptic  regime 
of  the  equilibrium  equations,  or  at  the  latest,  when  the  ellipticity  condition  fails.  Thus, 
localized  modes  of  deformation  at  the  boundary  or  at  the  interface  generally  occur  before 
the  onset  of  so-called  shear  banding  modes. 

3.3.  The  non-linear  case  :  some  results 

Although  the  complete  analysis  of  the  non-linear  problem  is  not  yet  available, 
some  results  can  be  given  for  the  possibility  of  emergence  of  deformation  modes 
involving  jumps  of  the  velocity  gradient  for  the  bi-linear  rate  constitutive  laws  (1). 

The  necessary  and  sufficient  conditions  for  the  onset  of  such  modes  inside  the 
body  have  been  given  by  Borre  &  Maier  who  extended  the  results  given  by  Rice,  and 
Rudnicki  &  Rice  for  so-called  continuous  and  discontinuous  localizations.  We  have 
amplified  these  results  by  seeking  necessary  and  sufficient  conditions  for  which  a 
discontinuity  surface  for  the  velocity  gradient  appears  at,  or  reaches  the  boundary  of  the 
solid.  These  conditions  are  given  below  for  the  constitutive  laws  (1)  with 
(3  :  z)  ®  (b  :  3) 

l~-  ~  h 

where  it  is  assumed  that  h  >  0,  and  3  is  strictly  positive  definite. 

At  a  point  P  of  the  boundary  T  where  only  surface  traction  rates  F  are  applied,  the 
necessary  and  sufficient  conditions  for  continuous  localization  [i.e.  the  material  is  in 
loading  (L  =  H)  on  each  side  of  the  singular  surface]  are 

•  •  • 

i)  there  exists  Sq  such  that  m  .  H  :  e0  =  F 

(2a)  ii  det  (n  .  K  .  n)  =  0 

iii)  (m  .  E  .  n)  .  (n  .  E  .  n)-1  .  (n  .  E  :  a)  =  m  .  E  :  a 
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At  a  point  P  of  the  boundary  F  where  only  surface  traction  rates  F  are  applied,  the 
necessary  and  sufficient  conditions  for  discontinuous  localization  [i.e.  the  material  is  in 
loading  (1  =  H)  on  one  side  and  in  unloading  (3  =  3)  on  the  other  side  of  the  singular 
surface]  are 

•  •  •  • 

i)  there  exists  So  such  that  m  .  13  :  £0  =  F  and  b  :  3  :  £0  >  0 

(2b)  ii)  det  (n  .  H  .  n)  <  0 

iii)  (m  .  3  .  n)  .  (n  .  3  .  nV1  .  (n  .  3  :  a)  =  m  .  3  :  z 

In  these  conditions,  m  denotes  the  unit  outward  normal  to  the  boundary  of  the  solid  in  P, 
whereas  n  denotes  the  unit  normal  to  the  singular  surface  in  P. 

At  a  point  Q  of  an  interface  I,  two  types  of  singular  surface  may  occur.  The 
singular  surface  can  either  stop  at,  or  cross  the  interface  (see  Fig. 2).  In  the  latter  case,  the 
singular  surface  can  meet  the  interface  at  different  angles  on  each  side.  Analogous 
conditions  to  conditions  (2)  can  be  exhibited. 

Remarks  : 

-  conditions  (2)  are  to  be  compared  with  the  corresponding  conditions  inside  the 
solid  given  by  Borre  &  Maier  and  stated  as 

(3)  det  (n  .  3  .  n)  <  0 

with  equality  corresponding  to  continuous,  and  inequality  corresponding  to  discontinuous 
localization  ; 

-  given  conditions  (2)  and  (3),  singular  surfaces  of  the  type  discussed  here 
generally  appear  First  inside  the  body  ; 

-  similar  conditions  to  conditions  (2)  can  be  exhibited  for  general  boundary 
conditions.  (Note  that  for  displacement  boundary  conditions,  relation  (3)  applies  both 
inside  the  body  and  at  the  boundary'.) ; 

-  conditions  (2>  and  (3)  are  a  priori  unrelated  to  the  boundary  complementing 
condition. 
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